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. Abstract 

The Fermat- Weber center of a planar body Q is a point in the plane from which the average 
distance to the points in Q is minimal. We first show that for any convex body Q in the 
■ plane, the average distance from the Fermat- Weber center of Q to the points of Q is larger 

than g • A(Q), where A(Q) is the diameter of Q. This proves a conjecture of Carmi, Har- 
Peled and Katz. From the other direction, we prove that the same average distance is at most 

2 ^13^ ' A(Q) < 0.3490 • A(Q). The new bound substantially improves the previous bound of 
^= • A(Q) w 0.3849 • A(Q) due to Abu-Affash and Katz, and brings us closer to the conjectured 

value of i • A(Q). We also confirm the upper bound conjecture for centrally symmetric planar 
convex bodies. 

^ ■ 1 Introduction 

The Fermat- Weber center of a measurable planar set Q with positive area is a point in the plane 
that minimizes the average distance to the points in Q. Such a point is the ideal location for a 
base station (e.g., fire station or a supply station) serving the region Q, assuming the region has 
uniform density. Given a measurable set Q with positive area and a point p in the plane, let hq(p) 
be the average distance between p and the points in Q, namely, 

J ggQ dist(p,g) dq 
= area(Q) ' 

where dist(p, q) = \pq\ is the Euclidean distance between p and q. Let FWq be the Fermat- Weber 
center of Q, and write fig = mhi{fjLq(p) : p G M 2 } = /jLq(FWq). 

Carmi, Har-Peled and Katz [3] showed that there exists a constant c > such that /ig > c-A(Q) 
holds for any convex body Q, where A(Q) denotes the diameter of Q. The convexity is necessary, 
since it is easy to construct nonconvex regions where the average distance from the Fermat- Weber 
center is arbitrarily small compared to the diameter. Of course the opposite inequality /j,q < c'-A(Q) 
holds for any body Q (convexity is not required), since we can trivially take d = 1. 
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Let ci denote the infimum, and C2 denote the supremum of fiq/A(Q) over all convex bodies Q 
in the plane. Carmi, Har-Peled and Katz [3] conjectured that c\ = g and 02 = \- Moreover, they 
conjectured that the supremum C2 is attained for a circular disk D, where /j,* D = | • A(D). They 
also proved that j < c\ < g. The inequality ci < g is given by an infinite sequence of rhombi, 
P £ , where one diagonal has some fixed length, say 2, and the other diagonal tends to zero; see 
Fig. [TJ By symmetry, the Fermat- Weber center of a rhombus is its center of symmetry, and one 
can verify that /Xp^/A(P e ) tends to |. The lower bound for c\ has been recently further improved 
by Abu-Affash and Katz from ^ to ^ pQ. Here we establish that c\ = g and thereby confirm the 
first of the two conjectures of Carmi, Har-Peled and Katz. 




Fi gure 1: A flat rhombus P e , with lim e _j.o Mp /A(-Pe) = g- 

Regarding the second conjecture, recently Abu-Affash and Katz proved that C2 < = 
0.3849.... Here we further improve this bound and bring it closer to the conjectured value of 
|. Finally, we also confirm the upper bound conjecture for centrally symmetric convex bodies Q. 

Our main results are summarized in the following two theorems: 

Theorem 1 For any convex body Q in the plane, we have fj,% > g • A(Q). 
Theorem 2 For any convex body Q in the plane, we have 

2(4 - y/S) 



< 
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A(Q) < 0.3490 • A(Q). 



Moreover, if Q is centrally symmetric, then fig < | • A(Q). 

Remarks. 1. The average distance from a point p in the plane can be defined analogously for 

finite point sets and for rectifiable curves. Observe that for a line segment / (a one-dimensional 

convex set), we would have /i}/A(I) = j. It might be interesting to note that while the thin rhombi 

mentioned above tend in the limit to a line segment, the value of the limit [f p / A(P £ ) equals |, not 
l 

A' 

2. In some applications, the cost of serving a location q from a facility at point p is dist K (p, q) for 
some exponent k > 1, rather than dist(p, q). We can define Hq(p) = (^f q( .Qdist K (p,q) dq^j /area(Q) 

and fi*Q = inf{/ig(p) : p G M 2 }, which is invariant under congruence. The ratio /ig*/A re (Q) is 
also invariant under similarity. The proof of Theorem [1] carries over for this variant and shows 
that /ig*/A K (Q) > ( K+ 2)2^ f° r an y convex body Q, and lim^o /ip*/2 K = ^, + ^ 2K . For the upper 
bound, the picture is not so clear: /ig/A(Q) is conjectured to be maximal for the circular disk, 



2 



however, there is a k > 1 such that /Xg*/A K (Q) cannot be maximal for the disk. In particular, if 
D is a disk of diameter 2 and it! is a convex body of diameter 2 whose smallest enclosing disk has 
diameter more than 2 (e.g., a regular or a Reuleaux triangle of diameter 2), then /xjp < //^*, for a 
sufficiently large k > 1. Let o be an arbitrary point in the plane, and let -D be centered at o. Then 
f q&D dist K (o, q) dq = J Q 27r r K ■ r dr d6 = and so lim^oo ^ < lim*^^ ^ = °- 0n tne 

other hand, for any region R' lying outside of D and for any n > 1, we have / 9g ^/ dist K (o, g) dg > 
area(-R') > 0. If /?' = \ D is the part of R lying outside D, then lim^oo /ug" > area(i?')/7r > 0. 

Related work. Fekete, Mitchell, and Weinbrecht [8] studied a continuous version of the problem 
for polygons with holes, where the distance between two points is measured by the L\ geodesic 
distance. A related question on Fermat- Weber centers in a discrete setting deals with stars and 
Steiner stars [5JE]- The reader can find more information on other variants of the Fermat- Weber 
problem in [IJ [IT] . 

2 Lower bound: proof of Theorem [L] 

In a nutshell the proof goes as follows. Given a convex body Q, we take its Steiner symmetrization 
with respect to a supporting line of a diameter segment cd, followed by another Steiner symmetriza- 
tion with respect to the perpendicular bisector of cd. The two Steiner symmetrizations preserve the 
area and the diameter, and do not increase the average distance from the corresponding Fermat- 
Weber centers. In the final step, we prove that the inequality holds for a convex body with two 
orthogonal symmetry axes. 

Steiner symmetrization with respect to an axis. Steiner symmetrization of a convex figure 
Q with respect to an axis (line) £ consists in replacing Q by a new figure S(Q,£) with symmetry 
axis £ by means of the following construction: Each chord of Q orthogonal to £ is displaced along 
its line to a new position where it is symmetric with respect to £, see \12\ pp. 64]. The resulting 
figure S(Q,£) is also convex, and obviously has the same area as Q. 

A body Q is x-monotone if the intersection of Q with every vertical line is either empty or 
is connected (that is, a point or a line segment). Every x-monotone body Q is bounded by the 
graphs of some functions / : [a, b] — > K and g : [a, b] — )■ 1R such that g(x) < fix) for all x G [a, b]. 
The Steiner symmetrization with respect to the x-axis £ x transforms Q into an x-monotone body 
S(Q,£ X ) bounded by the functions h(f(x) — g(x)) and ^(g(x) — f(x)) for x G [a,b]. As noted earlier, 
area(S , (Q, £ x )) = area(Q). The next two lemmas do not require the convexity of Q. 

Lemma 1 Let Q be an x-monotone body in the plane with a diameter parallel or orthogonal to the 
x-axis, then A(Q) = A(S(Q,£ X )). 

Proof. Let Q' = S(Q,£ X ). If Q has a diameter parallel to the x-axis, then the diameter is 
[(a,c), (b,c)], with a value c£l, g(a) = c = f(a) and g(b) = c = f(b). That is, A(Q) = b — a. 
In this case, the diameter of Q' is at least b — a, since both points (a,0) and (6,0) are in Q' . If 
Q has a diameter orthogonal to the x-axis, then the diameter is [(xo,/(xo)), (xo,g(xo))] for some 
xo 6 [a, b], and A(Q) = /(xo) — g(%o)- I n this case, the diameter of Q' is at least /(xo) — g{xo), 
since both points (xo, |(/(^o) — 9( x o))) and (xo, \{g{xo) — /(xo))) are in Q'. Therefore, we have 
A(Q') > A(Q). 

Let Ai and Ai be two points on the boundary of Q' such that A(Q') = dist(Ai, A2). Since Q' 
is symmetric to the x-axis, points A\ and A2 cannot both be on the upper (resp., lower) boundary 



3 



of Q' . Assume w.l.o.g. that A\ = (x\, ^(f(x\) — g(xi))) and A2 = (X2, \{g{x2) — /(X2))) for some 
a < Xi, X2 < b. 

A(Q') = dist^) = ^(x 2 -^+(^ +/(X2) - g(Xl) - g(S2 y. 

Now consider the following two point pairs in Q. The distance between B\ = {x\,f{x\)) and 
B2 = (x 2 ,y(x 2 )) is dist(.Bi, £> 2 ) = \J (x 2 — xi) 2 + (f(x\) — S r (x 2 )) 2 ■ Similarly, the distance between 
Ci = (x 1 ,g(x 1 )) and C 2 = (x 2 ,f(x 2 )) is dist(Ci,C 2 ) = V (^2 - xi) 2 + (g{xi) - /(x 2 )) 2 . Using the 
inequality between the arithmetic and quadratic means, we have 

^ f(x l ) + f{x2)-g{x l )-g{x2) \j i < U{x l )-g(x2)) 2 + {g(x 1 )-f{x2))\ 

This implies that d\st{A u A 2 ) < max(dist(Si, B 2 ), dist(Ci, C 2 )), and so A(Q') < A(Q). We con- 
clude that A(Q) = A(S(Q,£ X )). □ 



Lemma 2 If Q is an x-monotone body in the plane, then /j,q > H* S (q e y 
Proof. If (xq, yo) is the Fermat- Weber center of Q, then 



. _ faf 9 f $V(*-zo) 2 + (y-yo) 2 <iy 

^ Q area(Q) 



dx 



Observe that \/ i x — xo) 2 + {y — yo) 2 dy is the integral of the distances of the points in a line 
segment of length f(x) — g{x) from a point at distance \x — Xo| from the supporting line of the 
segment. This integral is minimal if the point is on the orthogonal bisector of the segment. That 
is, we have 



p/(x) _ -, , - - ^ . r m .l,. f^-^)Y dy 

lg(x) Jg(x) 

4(/(x)- 9 (x)) 
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±(s(x)-/(x)) 
Therefore, we conclude that 

^ area(Q) 

L ■%(.)-/(,)) V 0* " ^o) 2 + y 2 dy dx 



/ V( x - x o) 2 + y 2 dy- 



area(S(Q,x)) = "W^' ° )} " 



□ 
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Triangles. We next consider right triangles of a special kind, lying in the first quadrant, and 
show that the average distance from the origin to their points is larger than 3. 

Lemma 3 Let T a right triangle in the first quadrant based on the x-axis, with vertices (a, 0), 
(a, b), and (1,0), where < a < 1, and b > 0. Then ^t(o) > 3. 

Proof. We use the simple fact that the x-coordinate of a point is a lower bound to the distance 
from the origin. 

/a i (/o Ki-*)/(i-*) dy) dx /a i^(i-*)/(i-«) dy) dx 

m ° j " 6(l-o)/2 6(l-o)/2 



^/^(l-^dx 2 /x 2 



x 



1 



6(l-o)/2 (1-a) 2 V 2 3 

2 (2a 3 -3a 2 + l) 2 (1 - a)(l + a - 2a 2 ) 



(1-a) 2 6 (1-a) 2 6 

1 (1 + a - 2a 2 ) 1 
(1 - a) 3 -3' 

The last inequality in the chain follows from < a < 1. The inequality in the lemma is strict, since 
\J x 2 + y 2 > x for all points above the x-axis. □ 

Corollary 1 Let P be any rhombus. Then fi* p > | • A(P). 

Proof. Without loss of generality, we may assume that P is symmetric with respect to both the 
x-axis and the y-axis. Let us denote the vertices of P by (—1,0), (1,0), (0, —b), and (0,6), where 
b < 1. We have A(P) = 2. By symmetry, //p equals the average distance between the origin (0,0) 
and the points in one of the four congruent right triangles forming P. Consider the triangle T in 
the first quadrant. By Lemma [3] (with a = 0), we have \x P = Ht(°) > 3- Since A(P) = 2, we have 
H P > i • A(P), as desired. □ 

Lemma 4 Lei T be a triangle in the first quadrant with a vertical side on the line x = a, where 
< a < 1, and a third vertex at (1, 0). Then [Pr(o) > 

Proof. Refer to Fig. [2jii). Let U be a right triangle obtained from T by translating each vertical 
chord of T down until its lower endpoint is on the x-axis. Note that area(T) = area(£7). Observe also 
that the average distance from the origin decreases in this transformation, namely Ht{o) > fJ,jj(o). 
By LemmaEl we have ^u(o) > §, and so ht{o) > 3, as desired. □ 

We now have all necessary ingredients to prove Theorem [TJ 



Proof of Theorem [IJ Refer to Fig. [2j Let Q be a convex body in the plane, and let c,d G Q 
be two points at A(Q) distance apart. We may assume that c = (—1,0) and d = (1,0), by a 
similarity transformation if necessary, so that A(Q) = 2 (the ratio //q/A(Q) is invariant under 
similarities). Apply a Steiner symmetrization with respect to the x-axis, and then a second Steiner 
symmetrization with respect to the y-axis. The resulting body Q' = S(S(Q,£ x ),£ y ) is convex, and 
it is symmetric with respect to both coordinate axes. We have A(Q') = A(Q) = 2 by Lemma [TJ 
and in fact c, d G Q' . We also have jJL*Q, < /ig by Lemma [2j 
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Let Qi be the part of Q' lying in the first quadrant: Q\ = {{x, y) £ Q' : x, y > 0}. By symmetry, 
FWq> = o and we have /j,q, = /Uq'(o) = fj,Q 1 (o). Let 7 be the portion of the boundary of Q' lying 
in the first quadrant, between points b = (0, h), with < h < 1, and d = (1, 0). For any two points 
p,q £ 7 along 7, denote by ^{p, q) the portion of 7 between p and g. Let r be the intersection point 
of 7 and the vertical line x = 3. 

For a positive integer n, subdivide Qi into at most 2n + 2 pieces as follows. Choose n + 1 
points b = qi, q2 ■ ■ ■ , q n +i = ?~ along 7(6, r) such that is the intersection of 7 and the vertical line 
x = (i — l)/3n. Connect each of the n + 1 points to d by a straight line segment. These segments 
subdivide Q\ into n + 2 pieces: the right triangle To = Abod; a convex body Qo bounded by rd 
and 7(r, <i); and n curvilinear triangles Aqidqt+i for i = 1, 2, . . . , n. For simplicity, we assume that 
neither Qo, nor any of the curvilinear triangles are degenerate; otherwise they can be safely ignored 
(they do not contribute to the value of fJ-q,)- Subdivide each curvilinear triangle Aqidqi + i along the 
vertical line through qi + \ into a small curvilinear triangle Si on the left and a triangle Tj incident 
to point d on the right. The resulting subdivision has 2n + 2 pieces, under the nondegeneracy 
assumption. 

By Lemma O we have /j.t (o) > |. Observe that the difference /^t (o) — g does not depend on 
n, and let 5 = ht {o) — |. By Lemma HI we also have /^(o) > |, for each i = 1,2, ... ,n. Since 
every point in Qo is at distance at least | from the origin, we also have ^q (o) > |. 

For the n curvilinear triangles Si, i = 1,2, ... ,n, we use the trivial lower bound /J>Si(o) > 0. 
We now show that their total area s n = Y17=i area(Si) tends to if n goes to infinity. Recall 
that the y-coordinates of the points qi are at most 1, and their x-coordinates are at most |. This 
implies that the slope of every line qid, i = 1, 2, . . . , n + 1, is in the interval [—3/2,0]. Therefore, 
Si is contained in a right triangle bounded by a horizontal line through qi, a vertical line through 
5t+ij a nd the line qid. The area of this triangle is at most 3(3" " (| " ^7)) = l/(12n 2 ). That is, 
s n = J27=i area(Si) < l/(12n). In particular, s n < 5 ■ area(To) for a sufficiently large n. Then we 
can write 

f peQl dist(o,p)dp ^ Q0 ( o ) ■ area(Qo) + EHqMt^o) • area(T i ) 

VQii ) = 777^ > 7777 

area(C^ij area(C^i) 

g(area(Qi) — s n ) + 5 ■ area(T ) 1 2(5 • area(T ) 1 

" area(Qi) — 3 3-area(Qi) 3 
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This concludes the proof of Theorem [TJ □ 

Remark. A finite triangulation, followed by taking the limit suffices to prove the slightly weaker, 
non-strict inequality: fig > | ■ A(Q). 

3 Upper bounds: proof of Theorem [2] 

Let Q be a planar convex body and let D = A(Q). Let dQ denote the boundary of Q, and let 
int(Q) denote the interior of Q. Let ft, be the smallest disk enclosing Q, and let o and R be the 

center and respectively the radius of f2. Write a = 2 ^~[^ ■ By the convexity of Q, o £ Q, as 
observed in [TJ. Moreover, Abu-Affash and Katz [TJ have shown that the average distance from o 
to the points in Q satisfies 

fi Q (o) < -^j= ■ A(Q) < 0.3850 • A(Q). 

Here we further refine their analysis and derive a better upper bound on the average distance 
from o to the points in Q: 

Hq(o) < 2(4 ~ 3 ^ } • A(Q) < 0.3490 • A(Q). 

Since the average distance from the Fermat- Weber center of Q is not larger than that from o, we 
immediately get the same upper bound on C2- We need the next simple lemma established in [1]. 
Its proof follows from the definition of average distance. 

Lemma 5 |T]. Let Q\, Q2 be two (not necessarily convex) disjoint bodies in the plane, and p be 
a point in the plane. Then /U(Q lUQ2 )(p) < max(u Ql (p), u Q2 (p)). 

By induction, Lemma [5] yields: 

Lemma 6 Let Qi,Q2, ■ ■ ■ ,Qn be n (not necessarily convex) pairwise disjoint bodies in the plane, 
and p be a point in the plane. Then 

M(QiU...uQ n )(p) < max(/i Ql (p),...^ Qn (p)). 

We also need the following classical result of Jung [10] ; see also 

Theorem 3 (Jung [10]). Let S be a set of diameter A(S') in the plane. Then S is contained in a 
circle of radius 4= • A (5). 

By Theorem [3] we have 

-D<R<-^=D. (1) 

Observe that the average distance from the center of a circular sector of radius r and center 
angle a to the points in the sector is 

Jq ax 2 dx ar 3 /3 2r 



ax dx ar 2 /2 3 
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Proof of Theorem [2]. If o E dQ then Q is contained in a halfdisk G of Q, of the same diameter 
D, with o as the midpoint of this diameter. Then by ([2]), it follows that Hq{o) < | • D, as required. 

We can therefore assume that o G int(Q). Let e > be sufficiently small. For a large positive 
integer n, subdivide O into n congruent circular double sectors (wedges) Wi, . . . ,W n , symmetric 
about o (the center of Q), where each sector subtends an angle a = ir/n. Consider a double sector 
Wi = Ui U Vi, where Ui and Vi are circular sectors of fi. Let Xi C Ui, and Yi C be two minimal 
circular sectors centered at o and containing Ui n Q, and n Q, respectively: Ui n Q C Xj, and 
Vi H Q C y"j. Let Xj and be the radii of Xi and Yi, respectively. Let X[ C Xi, and 1"/ C 7j be 
two circular subsectors of radii (1 — e)xi and (1 — e)yi, respectively. Since o £ int(Q), we can select 
n = n(Q,e) large enough, so that for each 1 < i < n, the subsectors X[ and Y- are nonempty and 
entirely contained in Q. That is, for every i, we have 

X'iUYfC WiHQCXiUYi. (3) 

It is enough to show that for any double sector W = Wi, we have 

JimM(H/nQ)(o) < aD, 

since then, Lemma [6] (with Wi being the n pairwise disjoint regions) will imply that Hq{o) < aD, 

concluding the proof of Theorem [2j For simplicity, write x = Xi, and y = y^. Obviously the 

diameter of W n Q is at most D, hence x + y < D. We can assume w.l.o.g. that y < x, so by 

Theorem [3] we also have x < -7= • D. Hence so far, our constraints are: 

V 3 

< y < x < —= ■ D and x + y < D. (4) 
v3 

By the minimality of the disk Q, the convex body Q either contains three points q\,q2, Q3 on the 
boundary of f2 such that the triangle (71(72(73 contains the disk center o in the interior, or contains 
two points qi,q2 on the boundary of f2 such that the segment (71(72 goes through the disk center 
o. In the latter case, the segment (71(72 can be viewed as a degenerate triangle (71(7293 with two 
coinciding vertices q2 and (73. 

Let r be the radius of the largest disk centered at o that is contained in the convex body Q. 
Then r is at least the distance from o to the longest side of the triangle (71(72(73, say q\Q2- Since 
|(?i(?2| < D, \oq\\ = 1 052 1 = R, w e have 



r > ^R 2 -D 2 /4. 
Then the constraints in ((H) can be expanded to the following: 

y/R 2 - D 2 /A <y<x<R< D/V3 and x + y<D. (5) 
By the definition of average distance, we can write 

/ pG (W/nQ) dist (°'P) d P 



are&(W n Q) 

2 3 ^ _ 2 

a{l-e) 2 - (f + £) 3 (1 



< "2 3 ^ _ 2 3 _ _ x Z. I* f(\] 

.-/-■> o • /-. £ )2.( x 2 + y 2)- W 



Let 

= 3 " ^TF' and = 3 " (i- £ )M* 2 + ^ - (?) 
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Clearly for any feasible pair (x,y), we have 

lim/i(x,y,e) = f(x,y). 

It remains to maximize f(x,y) subject to the constraints in ([5]). We will show that under these 
constraints, 

/(X , B) <M_V!U. (8) 

Then _ 

lim/*(wriQ)(°) < )xmfi(x,y,e) = f{x,y) < — D, 

as required. 

We next verify the upper bound in ([8]). Throughout our analysis, we may assume that D is a 
fixed constant and x, y, and R are variable parameters. Substituting z = y/x in (|Z|), we have 

f(x, y) = g(x, z) = y • *** 2 - 
Then, taking the partial derivative of g(x, z) with respect to z, we have 



9 x _ 2x / 3z 2 1 + z 3 
q^9(x,z) - — ■ ^Tz 2 ~ (i + z 2) 



2x 3z 2 (l + z 2 ) - (1 + z i )2z 2x z(z 3 + 3z - 2) 



3 (1 + z 2 ) 2 3 (1 + z 2 ) 2 ' 

The cubic equation z 3 + 3z — 2 = has exactly one real root zq = (y/2 + l) 1 / 3 — (y/2 — l) 1 / 3 = 
0.596 .... Thus for a fixed x, the function g(x, z) is strictly decreasing for < z < zo and is strictly 
increasing for zq < z < 1. Therefore, by the upper bound that x + y < D and the lower bound that 
\J R 2 — D 2 /4 < r < y in (|5|), the function f(x,y) is maximized when y takes one of the following 
two extreme values: 



yi 



sjR 2 - D 2 /A and y 2 = D - x. 



By the inequality that x < R < D/V3 in ©, it follows that x + yi < R + ^JR? - D 2 /A < 
D j\f3 + -D/\/l2 < Z). Since rr + y 2 = -D, we have y\ < y 2 - 

Case 1. We first consider the easy case that y = y 2 - Then x + y = D, and we have 

2 x 3 + y 3 _ 2 (x + y) 3 - 3(x + y)xy _ 2 D 3 — 3Dxy 
^ X,V) ~ 3 ' x 2 + y 2 ~ 3 (x + y) 2 - 2xy ~ 3 ' D 2 — 2xy ' 

Substituting u; = xy, we tranform the function f(x,y) to a function h\(w): 



f(x,y) = hi(w) = - 



2 3Dw - D 



3 



3 2w — D 2 
The function hi(w) is decreasing in w because 

d^ ~ 3 ' \2w-D 2 ~ (2w - D 2 ) 2 



2 3D(2w - D 2 ) - 2(3Dw - D 3 ) _ 2 -D 3 

3 (2w - D 2 ) 2 ~ 3 ' (2w - D 2 ) 2 ~ 
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Thus /(x,y) is maximized when xy is minimized. With the sum x + y fixed at D, and under 
the constraint that x < R < D/yE in (J3J) , the product xy is minimized when x = -^D and 

y = ^1 — -^=^ -D. Thus we have 

/(x, y) < \ ■ ^ 2 + ^~^ 2 D = 2{A -f ] D = 0.3489 ... D. (9) 
~ 3 (X\ + (l- X\ 13 



Case 2. We next consider the casa^that y = y\. With y fixed, the function f(x,y) is maximized 
when x is as large as possible because 

y) = -- ( - X3 + Ij3 2x 

dx ' 3 \x 2 + y 2 (x 2 + y 2 ) 2 

2 3x 2 (x 2 + y 2 ) - (x 3 +y 3 )2x 

~ 3 (x 2 + y 2 ) 2 

2 x(x 3 + 3xy 2 - 2y 3 ) 



2 i n,2\2 



3 (x 2 + y 

2 x(y 3 + 3y 3 - 2y 3 ) 
" 3 ' (x 2 + y 2 ) 2 



Thus for y = \J R 2 — D 2 /4 and under the constraint that x < -R in ([5]), the function /(x, y) is 
maximized when x = R and y = \J R 2 — D 2 /4: = \J x 2 — D 2 /4. It follows that 



dx dy d^x 2 — D 2 /4 x 

dl = 1 and dR = dR = y/ xZ-Dy* = X/V ' 



Let h^{R) = f(R, \J R 2 — D 2 /4). We next show that h>2{R) is increasing in i?. Taking the 
derivative, we have 

d L , m 2 ^3x 2 ||+3y 2 ^ x 3 + y 3 ^ dx , o dy 



d^ 2 ^ ~ 3 ' x 2 + y 2 (x 2 +y 2 ) 2 \ 2x dR + 2y dR 

2 /3x 2 + 3y 2 (x/y) x 3 + y 3 , . . . . 

1 ■' ■' (2x + 2y(x/y)) 



3 \ x 2 + y 2 (x 2 + y 2 ) 2 

2 (3x 2 + 3xy) (x 2 + y 2 ) - (x 3 + y 3 ) (2x + 2x) 

3 (x 2 + y 2 ) 2 

2 (3x 4 + 3x 2 y 2 + 3x 3 y + 3xy 3 ) - (4x 4 + 4xy 3 ) 

3 (x 2 + y 2 ) 2 

2 (x 4 + 3x 2 y 2 + 3x 3 y + xy 3 ) - (2x 4 + 2xy 3 ) 

3 (x 2 + y 2 ) 2 

((l + y/x) 3 -2-2(y/x) 3 ). 



x 4 



3 (x 2 + y 2 ) 2 



1 This case, when x + y < D, has been mistakenly overlooked in the proof given in [6]. 
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Substituting z = y/x, we simplify the last factor (1 + y/x) 3 — 2 — 2(y/x) 3 in the resulting expression 
above to 

h 3 (z) = (1 + z) 3 - 2-2z 3 . 

To show that ^h2(R) > 0, it remains to show that h 3 (z) > 0. For < z < 1, the function h 3 (z) 
is increasing in z because 

^h 3 (z) = 3(1 + zf - 6z 2 = -3(1 - z) 2 + 6 > 6 - 3 > 0. 
Recall that x > y. If R < 3 ^ 4 13 ^ D, then we would easily have 



2 x 3 + y 3 2 x 3 2 2 „ 2(4 - >/3) 

- • — n o < - • — ^ = -X < -R < . 

3 x 2 + y 2 ~ 3 x 2 3 ~ 3 ~ 13 



which matches the upper bound in case 1. Now suppose that R > W-^ D. Then 



° /R < 3(1^7!) " d ^ = = V^WWi > Jl - (jjjj^gj) 2 /4 = 0.2955 . . . . 
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It follows that 

M:)>/,,| /I] =0.122(5... >U. 

hence 

±*.m > o. 

We have shown that the function h^(R) is increasing in R. Then, under the constraint that 
R < D/y/3 in fi2(R) is maximized when R = A^D. Correspondingly, f(x,y) is maximized 

when x = -k=D and y = —}=D. Thus 

V 3 V 12 

r i v 



y ) < | . ^) + U) n = ^r-D = 0.3464 . . . D, (10) 
3 fj^ + ( ^J) 5 



which is (slightly) smaller than the upper bound obtained in case 1. This proves the upper bound 



m 



Centrally symmetric body. Assume now that Q is centrally symmetric with respect to a point 
q. We repeat the same "double sector" argument. It is enough to observe that: (i) the center of Q 
coincides with q, that is, o = q; and (ii) x = y < | • D for any double sector W. By ([6]), the average 
distance calculation yields now 

2x 3 2x D 

M(VKnQ) (°J < 3(1 _ e)2 . x 2 " ^737)2 < ^737)2 » 

and by taking the limit when e tends to zero, we obtain 

A*Q(o) < -3 ' 

as required. The proof of Theorem [5] is now complete. □ 
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4 Applications 



1. Carmi, Har-Peled and Katz [3] showed that given a convex polygon Q with n vertices, and a 
parameter e > 0, one can compute an e-approximate Fermat- Weber center q £ Q in 0(n + 1/e 4 ) 
time such that fJ>Q(q) < (1 + £)[J*q- Abu-Affash and Katz [Tj gave a simple 0(n)-time algorithm 
for computing the center q of the smallest disk enclosing Q, and showed that q approximates the 
Fermat- Weber center of Q, with HQ(q) < ^j^q- O ur Theorems [[] and [2j combined with their 
analysis, improves the approximation ratio to about 2.09: 

12(4 - yg) , 

2. The value of the constant c\ (i.e., the infimum of /ig/A(Q) over all convex bodies Q in the plane) 
plays a key role in the following load balancing problem introduced by Aronov, Carmi and Katz [2]. 
We are given a convex body D and m points pi,P2, ■ ■ ■ ,Pm representing facilities in the interior 
of D. Subdivide D into m convex regions, Ri, R2, ■ ■ ■ , R m , of equal area such that Y^Li Mpi(-^i) 
is minimal. Here n Pi (Ri) is the cost associated with facility pi, which may be interpreted as the 
average travel time from the facility to any location in its designated region, each of which has the 
same area. One of the main results in [2] is a (8 + \/27r)-iactor approximation in the case that D 
is an n\ x n,2 rectangle for some integers ni,n2 £ N. This basic approximation bound is then used 
for several other cases, e.g., subdividing a convex fat domain D into m convex regions Ri. 

By substituting c\ = g (Theorem [JJ into the analysis in [2] , the upper bound for the approxi- 
mation ratio improves from 8 + v2~7r ~ 10.5067 to 7 + y2n ~ 9.5067. It can be further improved 
by optimizing another parameter used in their calculation. Let S be a unit square and let s £ S be 
an arbitrary point in the square. Aronov et al. [2] used the upper bound /j,s{s) < §\/2 ~ 0.9429. 
It is clear that max s gg Ms( s ) i s attained if s is a vertex of S. The average distance of S from such 
a vertex, say v, is fx s (v) = \ (\/2 + ln(l + y/2)) w 0.7652, and so /x s (s) < \ (y/2 + ln(l + \/2)) , 
for any s £ S. With these improvements, the upper bound on the approximation ratio becomes 
7 + ^ (y/2 + ln(l + V2)) w 9.0344. 

Acknowledgment. We are grateful to Alex Rand for stimulating discussions. 
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